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Quantum-Mechanical Position Operator and Localization in Extended Systems 
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We introduce a fundamental complex quantity, zl, which allows us to discriminate between a 
conducting and non-conducting thermodynamic phase in extended quantum systems. Its phase can 
be related to the expectation value of the position operator, while its modulus provides an appro- 
priate definition of a localization length. The expressions are valid for any fractional particle filling. 
As an illustration we use zl to characterize insulator to "superconducting" and Mott transitions in 
one-dimensional lattice models with infinite on-site Coulomb repulsion at quarter filling. 



Macroscopically, the fundamental property that dis- 
tinguishes an insulator from a conductor is that a steady 
current cannot flow at zero temperature Q| . Understand- 
ing the quantum nature of the insulating state and possi- 
ble phase transitions to a conducting state is a significant 
endeavor which started with the seminal work of Kohn 
H . Localization of the electronic wave function and the 
existence of a dielectric polarization density field arc two 
related features of the insulating ground state (GS). Ac- 
cordingly, one would expect that wisdom developed from 
the recent microscopic theory of polarization || could be 
exploited to establish a criterion for localization. 

In this Letter we introduce a complex quantity, 
zt,, which enables us to distinguish between conduct- 
ing (metal, superconductor) and non-conducting (band, 
Peierls, Anderson, and Mott insulators) states of mat- 
ter. Its phase corresponds to the GS expectation value 
of the position operator, intrinsically connected to macro- 
scopic polarization Q] , while its modulus provides an un- 
ambiguous definition of a localization length which can 
be described entirely in terms of the properties of the 
many-body ground state. In this way, zl generalizes the 
concept and quantity zn recently introduced by Resta 
and Sorella Q which, as we will show below, is valid 
only at integer particle filling. Furthermore, we show the 
usefulness of these ideas to characterize metal-insulator 
or superconducting-insulator quantum phase transitions 
in strongly correlated lattice models. 

To describe the intrinsic bulk properties of extended 
quantum systems it is almost mandatory to assume N in- 
teracting particles enclosed in a box subjected to periodic 
Born- von Karman (BvK) boundary conditions (BC), to 
reduce surface effects. On the other hand, the fact that 
its quantum state is defined in a non-simply connected 
manifold makes the center of mass of its wave function 
an ill-defined concept. In recent interesting papers 
the appropriate definition of the GS expectation value of 
the position operator (X) || and the localization length 
A H in one-dimensional (ID) systems satisfying BvK BC 
were discussed. The primary quantity was 



where \g) is the many-body GS of the system ((g\g) = 1), 



N 



is the sum of particle positions (X/N is 



the center of mass coordinate), and L is the number of 
unit cells in the system. (X) was defined from the phase 
of zjv, while A was specified from the modulus of zn, 
taking the thermodynamic limit (N, L — > 00, keeping the 
density no = N/L constant). The crux of these ideas 
is the expectation that, in the thermodynamic limit, zn 
can attain only two distinct values: zn vanishes for a 
conducting phase, while \zn\ — > 1 for an insulator. 

In microscopic calculations, translational invariance 
and BvK (or twisted) BC are convenient, and in most 
cases imperative, since the Hamiltonian can be diago- 
nalized separately in each Hilbert subspace of the many- 
body states which belong to the same irreducible rep- 
resentation of the translation group (or space group, in 
general). For systems which do not possess translational 
invariance, BvK (or twisted) BC are not an advantage 
over open BC, and for the latter, the traditional defini- 
tion of the expectation value of the position operator (i.e. 
the first moment of the modulus squared of the many- 
body wave function, |<?| 2 , or "dipole") is a well-defined 
concept. Thus, in the following we concentrate only on 
translational invariant systems. 

For these systems, excluding accidental degeneracies, 
I g) belongs to a well-defined irreducible representation of 
the space group. In particular, if the GS is non degen- 
erate (as assumed in Refs. |§,||), calling T the operator 
which translates one unit cell (taken as the unit of length) 
to the right, one has T \g) = e %K \g), where K is the total 
momentum of the GS. Degeneracies of states which differ 
in K or any other conserved quantum number do not af- 
fect the arguments given below because the Hamiltonian 
does not mix them. In other words, energy level cross- 
ings are allowed whenever the states involved belong to 
different symmetry sectors. 

In order for the matrix element in Eq. (Q) to be dif- 
ferent from zero, the operator e l ~ x , when decomposed 
into irreducible representations, should contain the triv- 
ial representation of the space group. This is not the 
case if no is not an integer. It is straightforward to check 
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(g\f tf e iS £* f*f\g) = e- l2lxn ^g\e^ x \g) , and the ma- 
trix clement should vanish unless uq is an integer. This 
restriction is too severe. For example, in ID a Mott tran- 
sition can take place for any filling, particularly if no is 
given by a simple fraction Jp| . For non- integer fillings, 
the definitions of Refs. |^,||| lead to an undefined (X) 
and infinite A for any system which is incorrect. 

The correct definition for arbitrary fillings no — n/l, 
where n/l is an irreducible fraction, is 

z L [n/l] = (g\e i ^ lx \g) . (2) 

Then, the GS expectation value of the position operator 
is defined (mod(L/7)) as 

(X [n/l]) = ^Imlnz^n/Z] , (3) 

and similarly to Ref. || one can introduce a localization 
length A. The first important point in these definitions 
is that the operator e l ~ lx is invariant under transla- 
tions and has no definite parity under space inversion. 
This ensures that the matrix element entering Eq. (|^) 
is not zero in finite systems, except for very particular 
cases in which hidden quantum numbers exist, like for 
free electrons, in which not only the total wave vector 
K, but also the one-particle wave vectors are conserved 
by the Hamiltonian, while the latter are shifted by the 
operator e l ^ lx |Q,|J. The second important remark is 
that, for non-integer no, when a family of Hamiltonians 
is introduced in which the one-particle wave vectors are 
shifted by a parameter a (Eq. (5) of Ref. g), there is 
a crossing of levels as a function of a between states of 
different if's. Crossings of these type were found previ- 
ously when calculating Berry's phases in systems which 
show the so-called anomalous flux quantization charac- 
teristic of superconductors Jt],^| , and are harmless, since 
the states involved have different quantum numbers and, 
therefore, are not mixed by the Hamiltonian. However, 
as a consequence of this crossing, Eq. (7) of Ref. jlj, 
stating that e %2 ^ x \g) = e llL \g) + 0(1/ L) (where 7l is a 
geometric phase) is incorrect, since both members have 
different if's for non- integer no. Finally, the behavior of 
zl in the thermodynamic limit provides us with a univer- 
sal criterion to distinguish between a conductor (metal or 
superconductor, as we will see below) and an insulator: 
zl vanishes in the first case, while \zl\ — > 1 in the second. 

In the general case (except for the above men- 
tioned particular cases with hidden symmetries) one has 
e i2 r lx \g) = e llL \g)+0(l/L), and the rest of the demon- 
stration leading to Eq. (||) can be done following Ref. 
Q with simple and straightforward changes. As shown 
elsewhere @, the Berry's phase jl — 2irl(X)/L is 

7 l = / d<z> (g K m\^g K ^)) , (4) 



where \gx (*&)) is the state obtained by adiabatic contin- 
uation of the GS \g), when a flux <E> is threaded through 
the ring |l(| (a ID system with BvK BC is topologically 
equivalent to a ring.) 

Using zl to distinguish between a conducting and a 
non-conducting system is more stringent than the crite- 
rion introduced by Kohn Kohn's criterion, univer- 
sally used in actual calculations, is based on the value of 
the charge stiffness D c (i.e. sensitivity of the curvature 
of the total energy to phase changes in the BC): in the 
thermodynamic limit D c is zero in the insulating phase, 
while it can attain any positive value in a conducting sys- 
tem. Instead, liniL^oo |zl| =0 for conducting systems, 
and Hin^^oo \z^\ = 1 for non-correlated insulators, or 
systems in which all particles are localized in non over- 
lapping sectors ||. To prove that this last limit can be 
achieved for any correlated insulator is rather involved, 
but one can prove, for a generic insulating state, that for 
large L, \zl\ — 1 — D c /n , and then lim^oo \z L \ = 1. 

To illustrate the main concepts described above, we 
consider an extended Hubbard model (EHM) with arbi- 
trary nearest-neighbor interaction V and infinite on-site 
repulsion U in ID, at quarter filling (no = 1/2). By a 
straightforward extension of the methods used in similar 
models Jill] , it can be shown that this is mapped onto 
a simple model of spinless interacting fermions jL2|, for 
which some analytical results can be obtained in the ther- 
modynamic limit. We explore the possibility of quantum 
phase transitions using our generalizations Zl and (X). 
Finally, we discuss two generalizations of this model. 

The Hamiltonian which describes the charge dynamics 
of the model in a ring of L = 2N sites with BvK BC 
(cj +L = ct) threaded by a flux $ is 

(5) 

The gauge transformation c| = e y */ L c| transforms 
H(Q) into a Hamiltonian H(&) in which the phase fac- 
tors disappear at the cost of introducing twisted BC 
cj +L = e 1 *^ except for fluxes $ = 2tt x integer p2| . 

For V 3> t, the GS can be obtained by perturba- 
tion theory. For t = it is two-fold degenerate be- 
tween the charge-density- wave (CDW) states |1) and |2), 
with the particles occupying every second site: |1) = 

n^4j+il°)> l 2 > = u f=i4j\°)- These states are mixed 
in N th order perturbation theory by 2N\ processes in 
which all particles hop either to the left or to the right in 
some order. For odd N Jig] , the effective matrix element 
is — r(N)(t/V) N cos(<I , /2), where r(N) is a real positive 
number with 1 < r(N)/2 < N\. Thus, the eigenstates are 
\g m = (|1) + \2))/V2 and |^($)) = (|1) - |2))/V2. 
The former (latter) is the GS for $ = ($ = 2tt). It 
is interesting to note what happens if $ is varied adia- 
batically from to 2n: the wave vector K = of the 
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GS of -ff ($) remains the same, while the wave vector K 
in the representation of H (<&) varies as K — K + <E>no 
in general R|. Since H(2tt) — H(0), after completing 
this cycle, |ffo(^)) has evolved from the GS of H(0) with 
-ft' = to its first excited state with K — tt. After another 
cycle |go(47r)) returns to the GS of H(Q) except for the 
geometrical phase Eq. (|4|). It is immediate to see that 
X\2) = Ef=r(2i)|2) = iV(AT+l)|2) ; X\l) = iV(JV-l)|l>. 
Then, e I "E LX |g ($)) = ±\g v ($>)), where the sign depends 
on the parity of (N + l)/2. This leads to a vanishing z/v, 
as defined by Eq. @), but to a zl[1/2] = 1. 

For attractive V and |V| ^> t, all particles group to- 
gether. For i = the degenerate GS is composed of 
the states \PS) 3 = fi\PS) (0 < j < N - 1), and 

\PS) = nf =1 c]|0). Then, j^PS), = (-l) N + 1 \PS) j 
leading to z L [l/2] = (-l) N+1 . 

For arbitrary values of V and t, the spinless model 
is equivalent (via a Jordan- Wigner transformation) to 
an XXZ model 2<y>a JaSfS" 1 with N spins up, and 
Jx — Jy — 2£, J z = V. This model was solved by the 
Bethe ansatz [[u]. From the solution one knows that 
there is a Mott transition at V = 2t and a transition 
to the phase segregated (PS) state at V = —2t. We 
have calculated the correlation exponent K p of the EHM 
in the conducting phase (— 2t < V < 2t), solving the 
Bethe-ansatz integral equations for the energy per site 
e(no) of the equivalent XXZ model H], and using the 
expression K p — 7Ta/ D c /(2& 2 e/dn 2 ) ), where the Drude 
weight (or charge stiffness) is given by the Bethe ansatz 
result D c = 2-nt sin/^/[8/i(7r — fi)] with /i = arccos(F/2t) 
p5[ . We obtain that K p > 1 (superconducting correla- 
tions dominate at large distances) for — 2t < V < V c , 
with V c = —\[2t within our accuracy (10~ 3 ). It is clear 
that this model provides a rich zero temperature phase 
diagram and an interesting laboratory to study localiza- 
tion and transitions to superconducting (SC) states. 

Unfortunately the Bethe-ansatz wave function is quite 
difficult to handle. Therefore, we obtained the GS \g) 
for finite systems with up to L = 16 sites by the Lanc- 
zos method. As usual, \g) was taken at the value of $ 
which minimizes the GS energy E g (<&). In the spinless 
model, this corresponds to $ = (BvK BC) for odd N 
and $ = 7r (antiperiodic BC) for even N, and in both 
cases K — (l2). This choice reduces the dependence 
of the GS energy and zl with size and also leads to a 
slightly more abrupt change in \zl\ near the transitions. 
The resulting \zl\ calculated with Eq. (||) is represented 
in Fig. ^ for various sizes, together with the correspond- 
ing numerical results for D c = (L/2)d 2 E g (<&) / d<S> 2 , and 
the above mentioned Bethe-ansatz result for D c in the 
thermodynamic limit (L — > oo). From the latter quan- 
tity, the conductor-insulator transitions at |V| = 2t are 
evident. However, the finite-size results always lead to a 
non-zero D c and except perhaps for the change in cur- 
vature near V ~ —2.5t there are no clear indications 



of any transition. The situation improves if the finite- 
size results are extrapolated using a polynomial in X/L. 
The extrapolation agrees very well with the exact result 
in the conducting phase, and suggests a transition near 

V = —2t. However, from the extrapolated D c no conclu- 
sions concerning the Mott transition at V = 2t can be 
drawn. Only the large size dependence near V ~ 3i is 
indicative of a charge gap. 

Near V = — 2t, as V increases, \zl\ decreases abruptly 
from values near one, to very small values, as expected 
for an insulator-conductor transition. For V ~ 2t, the 
change in \zl\ is rather smooth and the size dependence 
is small, but also the change of behavior becomes more 
abrupt with increasing size. The extrapolated values of 
\zl\ clearly show an abrupt transition near V — —2t, in- 
dicating that the system is a conductor (very small \zl\) 
for — 1.8f < V < 1.2t, and suggesting that it is an insula- 
tor for V ~ 3t. Notice that the convergence to the ther- 
modynamic limit for D c and \zl\ is, in absolute value, 
about the same. However, for any value of V, D c de- 
creases with system size while, in general, \zl\ decreases 
well inside the metallic phases (|V| < 1.5t) and increases 
well inside the insulating phases (|V| ~ 3i). Thus, the 
results for \zl\ are complementary to Kohn's D c criterion 
B, \zl\ providing a more useful measure. 

We have also studied the behavior of 7^ (Eq. (|])) near 
the transitions. Because of inversion symmetry, 7l = 
or 7r (mod(27r)). In contrast to previous cases we 
do not find a jump in 7^ at the transition from the SC to 
the PS regime (V = —2t), or at the opening of the charge 
gap (V = 2t). However, as in the case of attractive U 0, 
there is a jump from 7^ = n to 7l = as the dominant 
correlation functions at large distances change from the 
superconducting to the CDW ones. For 8 < L < 16 we 
obtain with four digits accuracy V c = —y2t in perfect 
agreement with the value obtained from the numerical 
solution of the Bethe-ansatz integral equations. Instead, 
the corresponding jump in the phase of Zl takes place for 
V c ~ — 1.55i. This suggests that using (X) = L^l/{2t:1) 
and Eq. @) one obtains a faster convergence for (X) to 
the thermodynamic limit than using Eq. (0) . In addition, 
since in the gapless metallic phases lim^^oo \zi\ = 0, it 
might be difficult to identify V c from the corresponding 
zero of zl for sufficiently large L. As in Ref. JtJ, we also 
find jumps in jl for other values of V, in particular for 

V = and V = t, without an obvious physical meaning. 
In the rest of this Letter, we discuss in more detail the 

expectation value of the position operator, Eq. (^), and 
illustrate its physical meaning using two other simple ex- 
amples. First, note that if \g) in Eq. (|J) is replaced by 
one of the CDW states |1) or \2) discussed above, the 
result is (X) = in both cases. This is reasonable, since 
|1) and 1 2) differ in a translation, and are thus physically 
equivalent in translationally invariant systems. This is, 
in general, the reason why (X) is defined mod(L/l). As- 
sume now that in the model Eq. rtq) one considers an 
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attractive V and adds a next-nearest-neighbor repulsion 
V such that V > \V\ > t. For even N = L/2, \g) con- 
sists of a sequence of two nearest-neighbor sites occupied, 
the next two empty, and so on. Specifically \g±- ! r/2{'&)) = 
(\A)±if\A))/V2, with \A) = n^- 1 4-, lC i., 2 jO>. In 



addition, X\A) = E^io ( 8 i + 3 )l^> = N ( N ~ l / 2 )\ A )- 
Then, e l ^ x \g ±7r/2 ) = -\g±n/2), and Eq. (§) gives 
(X) = L/4 = 2V/2 (mod(£/2)). This is consistent with 
the transfer of half of the particles of the system, either 
to the left or to the right, necessary to build the state 
(\A) ± it\A))/y/2 from (|1) ± »f|l))/V5. 

In the examples discussed above, the Hamiltonian had 
inversion symmetry. Under those circumstances it can be 
shown that the phase of Zl can only attain the values 
or 7r {zl is real). This statement implies that well-defined 
values of (X) different from and Lj (21) (mod(L/Z)) can 
only be achieved in the absence of inversion symmetry. 
A simple example in which this symmetry is explicitly 
broken is obtained replacing each site in Eq. (0) by an 
heteronuclear molecule. Specifically, in the spinless case 
consider a system described by the Hamiltonian 

3 3 

+E^' c M d i+^ c Wi4^)- ( 6 ) 



It is easy to see that for V, V',t' 3> t, and uq — 1/2 (one 
particle each two unit cells), \g) is a CDW with every sec- 
ond molecule singly occupied in the GS of the first two 
terms of Eq. (§). For this \g), Eq. (§), gives (X) = NP 
where P is the polarizability of each molecule and varies 
continuously as a function of t'/A [|. 

In conclusion, we have introduced a complex quantity 
Zl which is shown to display a qualitatively different be- 
havior for conductors than for insulators, thereby pro- 
viding a valuable criterion to distinguish between those 
states of matter. The criterion, which is valid for any 
fractional particle filling and for both, ordered and disor- 
dered systems, is complementary and in a sense sharper 
than that based on the value of D c [||. We have shown 
how to use Zl to characterize Mott and superconduct- 
ing transitions in models of strongly correlated quantum 
particles. Zl involves only the computation of the GS of 
the system, which for time-reversal symmetric Hamilto- 
nians is real valued, whereas D c needs at least two GS 
calculations for different BC and in one of them, time- 
reversal symmetry is broken. This seems innocuous for 
Lanczos studies, but it is not for stochastic approaches 
where the study of non-time reversal symmetric states 
adds an additional complication to the already infamous 
fermion sign problem. 
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FIG. 1. Drude weight D c (top) and modulus of zl defined 
by Eq. (Q) (bottom) as a function of V/t for the U — » +oo 
EHM at filling no = 1/2 for rings of various lengths L. Also 
shown are the corresponding quantities extrapolated to the 
thermodynamic limit using a cubic polynomial in 1/L, and 
the exact result for D c in the thermodynamic limit. The 
vertical lines separate the regions of phase segregation (PS), 
Luttinger liquid with K p > 1 (SC) or K p < 1 (C) and a new 
insulating phase at large V/t (I). 
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Figure 1: A. A. Aligia and G. Ortiz 



